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ABSTRACT 

We present the first results of a targeted survey carried out with the 305 m Arecibo telescope to 
detect Hi-line emission from galaxies at redshift z > 0.16. The targets, selected from the Sloan Digital 
Sky Survey database, are non-interacting disk galaxies in relatively isolated fields. We present here 
the Hi spectra and derived Hi parameters for ten objects detected in this pilot program. All are 
massive disk galaxies in the redshift interval 0.17—0.25 (i.e. 2—3 Gyr look-back time), with Hi masses 
A/hi = 3 — 8x10^° M© and high gas mass fractions (Hi- to - stellar mass ratios ~10— 30%). Our 
results demonstrate the efficacy of exploiting Arecibo's large collecting area to measure the Hi mass 
and rotational velocity of galaxies above redshift z = 0.2. In particular, this sample includes the 
highest redshift detections of Hi emission from individual galaxies made to date. Extension of this 
pilot program will allow us to study the Hi properties of field galaxies at cosmological distances, thus 
complementing ongoing radio synthesis observations of cluster samples at z^O.2. 

Subject headings: galaxies: evolution — galaxies: spiral — radio lines: galaxies — cosmology: obser- 
vations 



1. INTRODUCTION 

Because cold atomic hydrogen represents the reservoir 
for future star formation and the physics of the 21 cm 
transition is relatively simple, detection of the Hi line 
emission has the potential to contribute a key quanti- 
tative ingredient to models of how galaxies acquire and 
retain gas and convert it into stars. Of critical impor- 
tance is the tracking of the gas content over the red- 
shift interval between z ~ 1 and 0, i.e. over the period 
when the observed cosmic star formation r ate density 
has d eclined by an order of magnitude ( e.g., iLillv et all 
119961 iMadau et aLlll996t iBell et al.|[2005l . and references 
therein). However, our knowledge of the Hi 21 cm line 
emission properties of galaxies is currently limited to the 
very local universe, and only recently have technical im- 
provements in radio frequency receivers and back-ends 
made attempts to explore the Hi emission from galaxies 
at z > 0.1 in any practical sense. 

The first detection of Hi emission from a galaxy 
at cosmological distance (z = .1766) was made by 
IZwaan. van Dokkum. fc Verheijeiil ([2001) with the West- 
erbork Synthesis Radi o Teles cope (WSRT). More re- 
cently, IVerheiien et al.l (j2007| ) successfully carried out 
a pilot program with the WSRT to image z ^ 0.2 
galaxy clusters at 21 cm. They targeted the Abell 963 
(z = 0.206) and Abell 2192 (z = 0.188) clusters for 

^ Max-Planck-Institut fiir Astrophysik, D-85741 Garching, Ger- 
many; bcatinel@mpa-garching.mpg.de. 

2 Arecibo Observatory, HC3 Box 53995, Arecibo, PR 00612, 
USA. The Arecibo Observatory is part of the National Astronomy 
and Ionosphere Center, which is operated by Cornell University 
under a cooperative agreement with the National Science Founda- 
tion. 

^ Center for Radiophysics and Space Research and National As- 
tronomy and Ionosphere Center, Cornell University, Ithaca, NY 
14853, USA. 

Department of Physics, University of Washington, Seattle, WA 
98195, USA. 

Department of Astronomy, University of Washington, Seattle, 
WA 98195, USA. 



20x12 hours and 15x12 hours respectively, and detected 
Hi emission from 42 galax ies. The G i ant M etrewave Ra- 
dio Telescope was used by iLah et all ()2007t ) to constrain 
the Hi content of star forming galaxies at z = 0.24. They 
obtained an average Hi spectrum by co-adding the non- 
detections of 121 galaxies with known optical redshifts. 
However, their detection is at best marginal. 

Historically, Arecibo's huge collecting area has offered 
a critical advantage to Hi line studies of extragalactic 
objects, because sensitivity is of paramount importance 
to the detection of weak Hi signals. Recent improve- 
ments to the Arecibo instrumentation, both front-end 
and back-end, are allowing the practical exploration of 
the frequency range corresponding to Hi at z > 0.1. Most 
notably, the installation of the new, single-pixel L-band 
receiver ("L-wide") in February 2003 has given access 
to frequencies down to 1120 MHz (corresponding to red- 
shifted Hi emission at z — 0.27) with improved system 
noise and performance. Naturally, single-dish Hi obser- 
vations of intermediate redshift galaxies are challenged 
by a lack of angular resolution. As the density of galax- 
ies per unit solid angle increases with redshift, so also 
does the probability of finding multiple sources within 
the telescope beam. However, beam confusion can be 
broken kinematically, if high-quality optical redshifts are 
available for the sources within the beam. The avail- 
ability of imaging and spectroscopic data from the Sloan 
Digital Sky Survey (SDSS; York et al. 2000) has enabled 
the identification of sources suitable for long-integration 
Hi line observations with the Arecibo telescope. In this 
Letter, we report the first results of a pilot survey carried 
out at Arecibo to detect Hi emission from disk galaxies 
at z - 0.2. 

All the distance-dependent quantities in this work are 
calculated assuming Hq = 71 km s~^ Mpc"'^, flu = 0.27 
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2. SAMPLE, OBSERVATIONS AND DATA REDUCTION 

The initial goal of this survey was to explore the feasi- 
bility of detecting Hi line emission from galaxies beyond 
the redshift range of past Arecibo surveys, that is at red- 
shifts z >0.05. First observations successfully detected 
over 50 galaxies with 0.05 < z < 0.09, and the program 
was extended to higher redshifts with the availability of 
the new L-wide receiver. In the end, the frequency in- 
terval corresponding to 0.11 < z < 0.16 was excluded 
because of the particularly hostile Radio Frequency In- 
terference (RFI) environment. Here, we report results 
only for objects with z > 0.16; the full dataset will be 
presented elsewhere (B. Catinella et al., in preparation). 

The targets for Hi spectroscopy were extracted from 
the SDSS database according to the following criteria: 
(a) observable from Arecibo during night-time (to min- 
imize RFI); (b) redshift 0.16 < z < 0.27 and not corre- 
sponding to the frequency of known RFI within that in- 
terval. This redshift range corresponds to the 1120—1220 
MHz frequency interval and is dictated by the instrumen- 
tation (see below) ; (c) inclination i > 45° (for use in disk 
scaling relations) ; (d) presence of Ha emission with line 
width between 100 and 700 km s~^ and equivalent width 
larger than 5 A; (e) exponential disk profile; (f) non- 
interacting, undisturbed, disk appearance in the SDSS 
images. The SDSS images of the candidate targets were 
carefully inspected to exclude objects with interacting 
or peculiar appearance, and/or with Hi emission possi- 
bly contaminated by that of other, nearby objects lying 
within the Arecibo beam. We also discarded targets ly- 
ing in the vicinity of continuum sources that would cause 
ripples in the baselines. 

The Hi data for the z > 0.16 sample were collected dur- 
ing four observing runs, distributed between Fall 2003 
and Spring 2007. Including all the overheads, the to- 
tal time spent on this part of the survey was 280 hours. 
The observations were conducted in standard position- 
switching mode: each observation consisted of an on/off 
source pair, each integrated for 4 minutes (5 minutes in 
the last observing run), followed by the firing of a cal- 
ibration noise diode. The off observation of a "blank 
sky" region was performed by tracking exactly the same 
range of azimuth and zenith angle as was previously oc- 
cupied by the on source, as commonly done at Arecibo 
to optimize the bandpass subtraction. The on-source 
integration time per object varied between 80 minutes 
and 4.7 hours, primarily depending on the redshift and 
on the observed rotational velocity of the target (for the 
same Hi flux, slower rotators are easier to detect). The 
spectra were acquired using the L-band wide receiver, 
which operates in the frequency range 1120—1730 MHz. 
We also used a 1120-1220 MHz filter and a 750 MHz 
narrow band front-end filter to limit the impact of RFI 
on our observations. Although sensitivity requirements 
suggest long accumulations of data, RFI considerations 
dictate fast spectral dumping, because much RFI is of 
short duration. Thus spectra were recorded every sec- 
ond with 9-level sampling. Two correlator boards, each 
configured for 12.5 MHz bandwidth, one polarization, 
and 2048 channels per spectrum (yielding a velocity res- 
olution of 1.8 km s~^ at 1200 MHz before smoothing) 
were centered at or near the frequency corresponding to 
the SDSS redshift of the target. Minimizing the number 



of different board settings facilitates the identification of 
RFI. In order to better discriminate terrestrial sources 
from cosmic ones, we also chose not to apply any on- 
line Doppler correction for the motion of the Earth (the 
correction is applied during off-line processing) . 

The data reduction was performed using our own rou- 
tines based on the standard Arecibo data processing li- 
brary and will be described in full detail elsewhere. In 
summary, the data processing of each polarization and 
on-off pair includes the following steps: Hanning smooth- 
ing, bandpass subtraction, RFI excision, and flux cali- 
bration. A final spectrum for each of the two orthog- 
onal linear polarizations for each galaxy is constructed 
by the combination of good quality records (those with- 
out serious RFI or standing waves), weighting each by 
a factor 1/rms^, where rms is the root mean square 
noise measured in the signal-free portion of the spec- 
trum. The two polarizations are then averaged, yielding 
the final spectrum. After boxcar smoothing and baseline 
subtraction, the Hi-line profiles are ready for the mea- 
surement of redshift, rotational velocity and integrated 
Hi line flux. Our measurement technique is explained, 
e.g., in lCatinella. Havnes. fc Giovanellil (|2007l §2.2). 

3. RESULTS 

The total number of galaxies targeted by our z > 0.16 
pilot survey is 41 (this number does not include sev- 
eral other objects that were initially targeted but had 
to be abandoned after short integrations because of RFI 
or other problems) . Half of the targeted galaxies are sure 
detections with high-quality Hi profiles, a subset of which 
is presented in this Letter. Of the remaining half, only 
~l/3 are clearly non-detections; the remaining objects 
include marginal detections, uncertain profiles or other- 
wise problematic cases that are still being analyzed. 

Figure[T] (Plate 1) presents 10 of the best quality Hi de- 
tections of z~0.2 galaxies obtained at Arecibo, ordered 
by increasing redshift (as indicated next to each spec- 
trum). These galaxies are not necessarily the highest 
signal-to-noise detections, but have been chosen to span 
the entire 0.16 < z < 0.25 redshift interval covered by 
our pilot program. Each panel of Fig. [1] shows the SDSS 
image (1' size, which is equivalent to about 1/4 of the 
half power full width of the Arecibo beam at 1200 MHz) 
and the Hi spectrum of the targeted galaxy (the full 12.5 
MHz bandwidth is shown). The Hi-line profiles are box- 
car smoothed by 15 channels, corresponding to a veloc- 
ity resolution of 27 km s~^ at 1200 MHz, and baseline- 
subtracted. The spectra are labeled with two numbers, 
the Arecibo General Catalog (AGC; a private database 
maintained at Cornell University by MPH and RG) iden- 
tifier and the SDSS redshift of the galaxy. The frequency 
and heliocentric velocity (bottom and top x axis, respec- 
tively) corresponding to the SDSS redshift are indicated 
by a dashed line. The dotted line is at the redshift mea- 
sured from the Hi spectrum. The last three objects, with 
z > 0.21, represent the highest redshift detections of Hi- 
line emission from individual galaxies made to date. 

Tabic [T] lists SDSS and Hi parameters for the galaxies 
presented in Fig. [1] similarly ordered by increasing red- 
shift. The quantities are listed as follows: Columns (1) 
and (2) are the source identifiers in the AGC and SDSS 
databases. Column (3) is the SDSS redshift, ^sdss! the 
associated uncertainty on the last digit is shown within 
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brackets. Column (4) is the absolute /-band magnitude 
in the Cousins system, Mi (see § I3.1|) . Column (5) is the 
galaxy inclination to the line-of-sight in degrees, i, esti- 
mated from SDSS tabulated axis ratios in r-band. Col- 
umn (6) is the on-source integration time of the Arecibo 
observation, Ton, in minutes (only on scans that were ac- 
tually combined are counted). Column (7) is the redshift, 
z, measured from the Hi spectrum. The error on the cor- 
responding heliocentric velocity, cz, is half the error on 
the width, tabulated in the following column. Column 
(8) is the observed velocity width of the source line profile 
in km s~^, W50 (measured at the 50% of the peak lev- 
els), with its statistical uncertainty. Column (9) is the 
velocity width corrected for instrumental broadening and 
cosmological redshift only, W§q. No inclination or turbu- 
lent motion corrections are applied. Column (10) is the 
observed, integrated Hi-line flux density in Jy km s~^, 
F = J S dv, measured on the smoothed and baseline- 
subtracted spectrum. The statis tical errors are calcu- 
lated according to equation 2 of 'Spring ob et al.l ()2005l . 
hereafter S05). Column (11) is the signal-to- noise ratio 
of the Hi spectrum, S/N, estimated following ISaintoiigi 
()2007t l and adapted to the velocity resolution of our spec- 
tra. This definition of S /N accounts for the fact that for 
the same peak flux a broader spectrum has more signal. 
Column (12) is the base-10 logarithm of the Hi mass, 
Mhi, in solar units. The Hi masses of these galaxies 
vary between 3 and 8x10^" Mq and were computed via: 



Ml 



HI 



2.356 X 10^ 

TT~z 



Mpc 



JS dv 
Jy km s" 



(1) 



where dij{z) is the hmrinosity distance to the galaxy at 
redshift z. 

As both Figure [1] and Table [T] demonstrate, there is 
an excellent agreement between SDSS and Hi redshifts 
(the largest discrepancy is Az = 0.0004 for AGC 252580 
and AGC 224321. However, the difference drops to 
Az = 0.0001 for AGC 224321 when the SDSS emission 
line redshiflP is used instead. The emission line redshift of 
AGC 212941, 0.2238, is also in slightly better agreement 
with the Hi one). The importance of having accurate 
SDSS redshifts is exemphfied by AGC 242147, the high- 
est redshift detection, whose spectrum shows strong RFI 
at 1137 MHz. Because of the long integration required 
and because in the first hour or two of observation one 
can see the RFI but not the galaxy's line emission, it 
is essential to know in advance precisely where the Hi 
signal is expected to appear. 

In § [2] we noted that, at fixed redshift and Hi flux, 
galaxies that appear to be slower rotators (i.e., either in- 
trinsically slower or just viewed more face-on) are signifl- 
cantly easier to detect. This point is clearly illustrated by 
the two galaxies shown in the second row of Fig. [H AGC 
262033 and AGC 224321, which have among the largest 
and the smallest observed velocity widths in this sample, 
respectively. These are also the objects that required 
the longest and the shortest integrations, respectively, in 
spite of having almost the same redshift and Hi flux. 

The rotational velocities deprojected to edge-on view 
(i.e., 0.5 Wgo/sin i) vary between ~180 and 320 km s~^ 
for this sample. These values are consistent with those 

<^ Not available for AGC 252580. 



typical of local disk galaxies with similar rest-frame I- 
band luminositie s, as shown for a large set of Ha ro- 
t at ion curves by ICatinella. Giovanelli. fc HavnesI (|2006l 
Fig. 1), also taking into account that optical velocity 
widths tend to be systematically smaller than 21 cm 
ones (a bias understandable in terms of rotation curve 
sha pes and relative extent of Hi vs. Ha emission; see, 
e.g.. ICatinella. Havnes. fc Giovanellill2007D . 

We computed gas mass fractions, Mni/Mstar, for the 
galaxies in our sample. Stellar masses, Mgtar, were es- 
timated from SDSS i-band mo del magnitudes and g — i 
colors, using iBell et el\ (|2003f l stellar mass-to-light vs. 
color correlations'". We corrected the S PSS magnitudes 
for Ga lactic attenuation first, and applied [Springob et al.l 
(|2007l hereafter S07) K- and internal extinction correc- 
tions (the extinction in SDSS g, i bands was obtained 
from that in S07 Cousins / assuming a standard attenu- 
ation curve of the form t\ cx A~° ''). The objects detected 
at Arecibo are very massive, gas-rich galaxies, with gas 
mass fractions Mni/Mstar of ~10— 30%, refiective of the 
deliberate selection of galaxies on the basis of their po- 
tentially high Hi content. 

3.1. Comparison with z = field galaxies 

Since z ^ 0.2 corresponds to a non-negligible fraction 
of the age of the universe, and in particular, the epoch 
over which the star formation density has rapidly de- 
clined, it is interesting to investigate if this small sam- 
ple is already suggestive of possible evolutionary effects. 
This exercise requires a representative comparison data 
set of field galaxies at z = which we have extracted 
from the Cornell Hi digital archive (S05). In addition to 
being a high-quality, homogeneous compilation of Hi pa- 
rameters for a large number of galaxies, a subset of the 
objects in the Hi archive has matching /-band photom- 
etry and grouping information from the SFIH — h catalog 
(S07). The cross-match of the S05 Hi archive (8850 ob- 
jects) with the S07 photometry (Table 2: 4053 objects, 
2713 of which not listed as group members), after rejec- 
tion of galaxies with recessional velocities cz > 10000 
km s~^, includes 1081 objects. The redshift cut was ap- 
plied in order to obtain a clean z = sample, but the 
inclusion of the 124 galaxies with cz > 10000 km s~^ 
(the largest recessional velocity is 17897 km s~^) would 
not change any of our conclusions. Since color informa- 
tion is not available for the SFI++ galaxies, we could not 
estimate stellar masses. Thus we compared Hi mass to /- 
band luminosity ratios. Mm/ Li, for the two samples. We 
used standard transformations* to convert SDSS model 
r magnitudes into Cousins / for our sample, and rescaled 
the SFI-I--I- quantities to our adopted Hubble constant. 
The resulting Mhi/Li ratios for our data set (large cir- 
cles) and for the z — reference sample (small symbols) 
are presented in Fig. ^ Dotted lines indicate lines of 
constant Hi mass, A/hi= 2, 4, 8 x 10^° M0 (bottom to 
top, respectively). The reference sample has a very large 
scatter, and contains only a small number of galaxies 
with rest-frame luminosities comparable with those of 



^ As IBell et al.l I I2003I) note, their choice of stellar initial mass 
function yields "the maximum possible stellar M/L ratio". Thus, 
at given luminosity, we might overestimate the stellar masses and 
therefore underestimate the gas mass fracti ons. 

* See http ://www.sds s. org/dr6/algorithms/sdssUBVRnYansform.html\ 
(Lupton 2UU5 equations J! 
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our sample. This limitation makes it very difficult to 
estimate where the center of the distribution might be, 
i.e. what is the "typical" Hi mass of a local field galaxy 
with Li~10^"'^ L0. The objects in our sample are, on 
average, very gas-rich, and lie on or slightly above the 
upper envelope of the local galaxy distribution. Natu- 
rally, this selection effect arises because we detect only 
the gas-richest galaxies at z = 0.2. Admittedly, galaxies 
with such large Mfji/ii ratios at these luminosities do 
appear to be rare in the local universe. However, since 
our survey probes a much larger volume compared to 
the z = sample, it is perhaps not surprising to detect 
several of these locally rare objects. 

Understanding and quantifying evolutionary effects re- 
quires the use of complete samples with well-defined se- 
lection criteria. Clearly, this is not the case for our cur- 
rent survey, which was in part a proof-of-feasibility ex- 
periment. Unfortunately, the z = reference sample is 
not complete in any sense, either. In particular, it is cur- 
rently very difficult to find local Hi samples with large 
numbers of galaxies that are as massive and/or as Hi- 
rich as our detections. Thus, we cannot make any strong 
statements about the evolution of the Hi mass-to-light 
ratio based on the results shown in Fig. [21 

4. CONCLUSIONS 

In this Letter, we have shown that carefully-planned 
Arecibo observations can deliver high-quality global 
Hi profiles for galaxies at redshifts up to z = 0.25 
with integration times of a few hours per object. The 
detected objects are massive disk galaxies with high 
gas mass fractions ('^10—30%), and are apparently 
rare in the local universe. At present, the small size of 
our sample, selection effects and the lack of adequate 
comparison samples at z = prevent us from drawing 
any firm conclusions about the evolution of the Hi 
content of galaxies between redshifts 0.25 and 0. Larger 
and well-defined samples of galaxies at z ^ 0.2 and 
higher are needed in order to make progress in this 
area. Equally importantly, fair comparison samples at 
z = are also needed. The latter might be provided 
by two large surveys that are ongoing at Arecibo. The 
first on e is the Arecibo Legacy Fast ALFA (ALFALFA) 
survey (jGiovanelh et al.ll2005l l2007| ). a bhnd Hi survey 
of the extragalactic sky visible from Arecibo, which 



will make a complete census of Hi-rich galaxies with 
z < 0.060 and has already detected a significa nt pop- 
ulatio n of objects with Hi masses > 10^° Mq (|Havnesl 
l2008h . The seco nd one is the GALEX Arecibo SPSS 
Surve y (GASS; ICatinella. S chimino vich. fc KauffmannI 
l2008f l. a recently started program designed to measure 
the Hi content of ^1000 massive galaxies, randomly 
selected from the intersection of the SDSS spectroscopic 
survey, GALEX and ALFALFA footprints based on 
their redshift (0.025 < z < 0.05) and stellar mass 
(10.0 < Log [Afstar/M©] < 11.5) only. ALFALFA and 
GASS will deliver unbiased samples of local Hi-rich and 
massive galaxies, respectively. In combination with the 
future extension of our pilot intermediate redshift field 
galaxy program and with complementary Hi synthesis 
mapping programs of similar galaxies in crowded cluster 
and group fields, it is clearly now feasible to explore the 
evolution of the Hi content of disk galaxies over the last 
several Gyr for a statistically complete sample. 
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Fig. 1. — SDSS postage stamp images and Hi spectra of the 10 Arecibo detections discussed in this work. Each SDSS image is 1 
arcmin square. The Hi spectra are cahbrated, smoothed and baseline-subtracted. Each spectrum is labeled with its AGC identifier and 
SDSS redshift. The central heliocentric velocity (frequency) corresponding to the SDSS redshift and that measured from the spectrum are 
indicated with a dashed and a dotted line, respectively. The rms noise per channel at the final velocity resolution (27 km s~^ at 1200 
MHz), measured in the RFI-free portions of the baseline, varies between 0.1 and 0.2 mJy for these spectra. 
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Log Li (Lq) 

Fig. 2. — Hi mass to /-band luminosity ratio as a function of /-band luminosity for the Arecibo Hi detections (large circles) and for a 
comparison field sample at 2 = (small symbols). Dotted lines correspond to constant Hi mass, A/hi= 2,4,8 X lO'^" Mq (bottom to top, 
respectively) . 
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SDSS AND Hi Parameters of the Arecibo Detections 
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J120948. 14+100822.5 


0.1879(2) 


-23.56 


46 


80 


0.1875 


417±14 


328 


0.26+0.05 


8.9 


10.63 


101750 


J003610.70+142246.4 


0.1909(1) 


-23.67 


62 


156 


0.1909 


542±21 


432 


0.19+0.03 


8.2 


10.51 


242091 


J140522. 72+052814.6 


0.1954(2) 


-24.00 


63 


156 


0.1953 


502±04 


397 


0.25+0.04 


10.3 


10.64 


232055 


J134211.36+053211.5 


0.2031(1) 


-23.45 


57 


92 


0.2030 


464±10 


362 


0.26+0.05 


8.4 


10.69 


262029 


J160938. 00+312958.5 


0.2188(2) 


-24.20 


53 


136 


0.2186 


642±56 


503 


0.34+0.05 


9.0 


10.88 


212941 


J111645. 15+054210.0 


0.2239(1) 


-23.84 


54 


180 


0.2236 


383±39 


289 


0.18+0.04 


7.0 


10.62 


242147 


J142735. 69+033434.2 


0.2454(1) 


-24.11 


51 


176 


0.2454 


499±11 


377 


0.28+0.05 


8.0 


10.90 



